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ABSTRACT 

We study the reliability of the reconstruction method which uses a modelling 
of the redshift distortions of the two-point correlation function to estimate the 
pairwise peculiar velocity dispersion of galaxies. In particular, the dependence 
of this quantity on different models for the infall velocity is examined for the Las 
Campanas Redshift Survey. We make extensive use of numerical simulations 
and of mock catalogs derived from them to discuss the effect of a self-similar 
infall model, of zero infall, and of the real infall taken from the simulation. The 
implications for two recent discrepant determinations of the pairwise velocity 
dispersion for this survey are discussed. 



Subject headings: Galaxies: formation - Cosmology: observations - dark matter 
- large-scale structure of Universe 



1. Introduction 



The Pairwise peculiar Velocity Dispersion (PVD) of galaxies is in principle a well- 
defined statistical quantity which can give interesting information on the cosmic matter 
distribution in addition to the two-point correlation function. The peculiar velocities of 
galaxies are determined by the action of the local gravitational fields, and thus they directly 
mirror the gravitational potentials caused by dark and luminous matter. The PVD is 
measured by modelling the distortions in the observed redshift-space correlation function 
£,z(r p ,Tr) which is in general not just a function of the distance s 2 = r 2 + 7r 2 , but depends 
anisotropically on the separations of a galaxy pair perpendicular to (r p ) and along (ir) the 
line of sight. This is the information that can be obtained from a redshift survey. The 
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basic step in modelling is to write £ z ( folding integral of the real-space correlation 

function £(r) and the distribution function f(vu) of the relative velocity t>i 2 of galaxy pairs 
along the line of sight 



1 + Ur P , tt) = y /(«ia) [l + eC^p + ^-W^o) 2 )] ^12, (1) 

(see, e.g., Davis & Peebles |1 983| ). The real space correlation function £(r) must be estimated 
from the redshift catalog through the relation 

poo poo , 

w ( r p)=J Ur P ,*)dn = J o ^r 2 p + y 2 )dy (2) 
where w(r p ) is the so-called 'projected' two-point correlation function. In most previous 



works a power-law form is assumed for £(r). Based on observational (Davis & Peebles |1983 
Fisher et al |1994| ) and theoretical considerations (Diaferio & Geller |1996| ; Sheth |1996| ; Seto 
& Yokoyama |1998|) an exponential form is usually adopted for f{y i 2 ): 



f{v 12 ) = —f= — exp \v 12 - v 12 \ , (3) 

where is the mean and <Ti2 is the dispersion of the 1-D pairwise peculiar velocities along 
the line of sight. It is worth pointing out that every step in the above modelling [eqs. (H]-|3|)] 
is only an approximation and the infall vu is unknown. As demonstrated by Jing, Mo, & 
Borner (|1998| ; JMB) (see also below) with mock samples of the Las Campanas Redshift 
Survey (LCRS; Shectman et al. |1996|) , however, the above procedure can give an accurate 
estimate of o\i (within 20% accuracy) if the infall is known. 

The distribution function j[v\2) is determined by its first and second moment: the 
infall (^12) and the dispersion (Jyi[t\ Both distort the two-point correlation function but in 
opposite ways. The infall velocity must also be modelled in some detail to allow a precise 
measurement of the dispersion Oyiir)- The situation seems somewhat complex: v^i in the 
real Universe is not known at the present. One might think that on small scales V12 is 
negligible, but this is true only for very small scales indeed. As has been shown by Mo, 
Jing, & Borner ( 1997| ) the function rises quite sharply around 1 h^ 1 Mpc, reaching twice 



the Hubble velocity just beyond 1 h 1 Mpc. Therefore it is necessary to model the infall 
carefully when measuring a^ir). 

In JMB we have determined the PVD for the LCRS using this reconstruction method. 
Recently an attempt to measure the PVD for the same survey directly from a Fourier 
deconvolution of the anisotropies of the redshift space two-point correlation function 
(Landy et al |1998| ) has resulted in a much lower value for the PVD (363 ± 44kms~ 1 vs 



570 ± 80kms 1 ). We think it is important to identify the causes for the discrepancy, since 
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the PVD is a very powerful test for the theories of the structure formation. We shall show 
that the different infall models used in the two studies can explain the discrepancy. But 
in addition to this immediate aspect there is the general question of how reliable these 
methods to measure the PVD really are. In this paper we want to address this question. 



2. N-Body Simulation and mock samples 

The true PVD can be easily determined from the three-dimensional velocities of 
particles in numerical simulations. Writing the three-dimensional velocity difference of 
particle pairs at points x and x + r, i.e. at separation r, as 

Vi2{r) = v(x) — v(x + r) (4) 

the true PVD is defined as 

a 12 (r) = ((v l2 (f) - (^ 12 (r-))) 2 /3) 1/2 (5) 

where (• ■ •) denotes the average over all pairs at separation r. In JMB we have considered 
three spatially flat cosmological models. Here we make use of one simulation of the model 
with density parameter Q = 0.2, cosmological constant Ao = 0.8, shape parameter T = 0.2 
and normalization ag = 1. This N-body simulation was generated with a P 3 M code (Jing 
and Fang 1994|) with 128 3 particles. Twenty mock samples without the fiber collisions are 



generated from the simulation to mimic the LCRS. We use these mock samples to determine 
the PVD in the way outlined above, and compare the results to the true PVD obtained 
from the simulations. Further details about the simulation, the mock samples, and our 
statistical method were given in JMB. 



3. Dependence on Infall-models 

In JMB we have determined the PVD for the LCRS and found a best estimate of 

cri 2 (l/i _1 Mpc) = 570 ±80 km s" 1 (6) 

at a separation of 1 /i~ 1 Mpc. Figure |T] gives an indication of the reliability of the result: 
accepting the exponential shape of the distribution function f(vu) as a reasonable ansatz, 
there remains the mean value vu(r) to be modelled. The filled squares in Fig. [IJ represent 
the results from the modelling process, when a self-similar infall model is used for v^: 

1 + (r/rj 2 
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and y is the radial separation in real space. This form for v& nas been widely used in 
previous work. It also is a good approximation to the real infall pattern in some CDM 
models. The circles in Fig. [I] denote the PVD reconstructed when the infall is set to zero. 
In Fig. 0, we have plotted the true PVD read off directly from the simulation as the solid 
line. Circles are the result from the mock catalogs with zero infall; triangles depict the 
result from the modelling process when the real infall from the simulation is used, and 
filled squares are again the reconstruction for a self-similar infall. The two different infall 
models give very similar results on scales r p < 10 h Mpc. We may conclude that the 
reconstruction of the PVD does not depend very sensitively on the model for v±^. The 
PVD reconstructed from the redshift distortions agrees qualitatively with the true value. 
There are, however, differences of some significance, even if the real infall pattern is used, 
which reflects the approximate features of the modelling [eqs. ([T]-[3])]. The model, where the 
infall is completely neglected, in contrast does not even qualitatively correspond to the true 
value. Although the infall velocity becomes negligible on small scales it still has a strong 
influence on the PVD: at 1 h~ 1 Mpc we find a reduction from 570 km s -1 to 400kms _1 , if 
Vx2 = 0, and at larger scales the no-infall model lets the PVD go down quite rapidly. The 
same behaviour of rapid drop of the no-infall models can be seen in the simulation results 
(Fig. |2]). (The simulation results are higher, because we have not yet applied the bias 
necessary to achieve agreement with the observation). It is probably due to the shape of 
the infall velocity (Mo et al. |1997| ) around 1 /i _1 Mpc with a steep rise, and a maximum at a 
few /i _1 Mpc for all CDM models considered in Mo et al. Q1997 ). The Fourier deconvolution 
method as applied by Landy et al. (|1998| ) appears simpler, because it does not seem to 
have to model the infall W^. In fact, however, they use the model ujj — which can, as we 
have seen from our simulations, lead to an underestimate of the true value. To obtain the 
true PVD therefore also the Fourier deconvolution method needs to model the infall, and 
thus it meets the same difficulties as the usual approach. 

The influence on the PVD estimate at ~ 1 /i _1 Mpc of the infall models we found 
here is in qualitative agreement with many previous works, e.g, Davis & Peebles (1983) 
and Marzke et al. (|1995 ). Quantitatively the influence may depend on the sample used, 
which is particularly true for small samples, since the infall effect depends on which regions 
(clusters or fields) the sample has surveyed. Certainly the infall effect should be universal 
for fair samples, but it is not known that any observational sample available to date could 
be considered "fair" for the infall effect. Therefore it is necessary to quantify this effect 
individually for each sample, as we have done for the LCRS here. 



- 5 - 



4. Discussion 



The modelling of the redshift distortions of the two-point correlation function gives 
a reasonable estimate — certainly within 20% of the true PVD, despite the complex 
reconstruction method involved. The differences are due to the fact that the form of f{v\2) 
is only approximately an exponential, and that the PVD estimated from the redshift 
distortions is some kind of average of the true PVD along the line of sight. Since the 
true PVD depends on the separation of galaxy pairs in real space, these two quantities 
are different by definition. Whether we use the self-similar or the true infall model has 
little effect on the results. It is very important , however, to use both the first and second 
moments of the velocity distribution function in the modelling process, since they lead 
to distortions of the redshift space correlations in opposing directions. Thus, for instance 
setting v^2 = leads to drastic changes. a X2 at 1 /i -1 Mpc drops from 570 to about 400 
km/s, and it becomes very small even for > 1 /i _1 Mpc. We suspect that this behaviour is 
responsible for the result of a recent work (Landy et al |199S| ), where a 12 from the LCRS is 



estimated to be 363 ± 44kms _1 at scale ~ 1 /i _1 Mpc. Their reconstruction method makes 
use of a Fourier deconvolution of £ z (r p , tc) , but the effects of infall are neglected. We can 
reproduce their result, if we set = , but from our comparison between simulations and 
mock catalogues we can draw the conclusion that this approach may have underestimated 
the true PVD and that the analysis incorporating a reasonable infall, like that of JMB, 
gives a much more reliable result. 

Considering the fact that the two works have used rather different methods to measure 
the PVD, we would stress that one should remain open minded with respect to both values 
even though we have quantitatively explained the discrepancy with the different infall 
models. However, the PVD is a very important quantity in cosmology. Since the procedure 
of JMB has been extensively tested with mock samples, it is important and necessary to 
make a similar test to the procedure of Landy et al. (|1998| ). 
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Fig. 1. — The PVD of the Las Campanas Redshift Survey: filled squares for the self-similar 
infall and circles for no infall. Error bars are la deviations given by bootstrap resampling. 
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Fig. 2. — The PVD of 20 mock samples without fiber collisions. Three infall models are 
adopted for Vu(r): the self-similar infall model (filled squares), the zero infall model (circles), 
the infall derived directly from the simulations (open triangles). The true pairwise velocity 
dispersion given by the 3-dimensional velocities in the simulations is shown as the solid line. 
Error bars are the (la) standard deviations of the mean from the mock samples. 



